Development of Knife-Edge Ridges on Ion-Bombarded Surfaces 
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We demonstrate in both laboratory and numerical experiments that ion bombardment of a modestly sloped 
surface can create knife-edge like ridges with extremely high slopes. Small pre-fabricated pits expand under 
ion bombardment, and the collision of two such pits creates knife-edge ridges. Both laboratory and numerical 
experiments show that the pit propagation speed and the precise shape of the knife edge ridges are universal, 
independent of initial conditions, as has been predicted theoretically. These observations suggest a novel 
method of fabrication in which a surface is pre-patterned so that it dynamically evolves to a desired target 
pattern made of knife-edge ridges. 



The efficient fabrication of ever-smaller structures is 
one of the major challenges of 21st century science and 
engineering. Ion bombardment has emerged as a promis- 
ing candidate to create patterns on surfaces 14 . One 
method uses a focused ion beam to micro-machine sharp 
features directljEHI] _ this allows for detailed control of 
the shape of the features but, as a serial writing process, 
is too time-consuming to pattern large areas. Another 
method is to bombard a surface uniformly, which can 
excite linear instabilities that grow into patterns such as 
quantum dots^E] This is less costly, but can only create 
structures as small as the smallest linearly unstable wave- 
length, with steepnesses limited due to saturation of the 
linear modes. To overcome these limitations one would 
like to create steep, sharp structures spontaneously, by 
exploiting the dynamical processes that underlie surface 
evolution under ion bombardmenlP^l. 

Recently we proposed a scenario for creating very 
sharp features on ion bombarded surfaces, by starting 
with a surface that is pre-patterned to have modest 
slopes on the macroscale^. Under ion bombardment, 
our theoretical calculations demonstrated that if the ini- 
tial slope is in the right range, the structures would spon- 
taneously evolve to knife-edge-like ridges, with extremely 
high slopes, and high radius of curvature. Both the final 
slope and radius of curvature are independent of initial 
conditions, and depend only on the shape of the curve 
describing sputter yield (atoms out per incident ion) vs. 
incidence angle. Here we demonstrate the formation of 
knife edge ridges in experiments. Our experiments show 
that uniformly irradiated small pits expand outward, de- 
veloping steep sides with uniform slopes. When two pits 
collide, the front evolves to a sharp, knife-edge-like struc- 
ture with features on a scale much smaller than any con- 
tained in the initial conditions. Numerical simulations 
of the classical macroscopic equations show remarkably 
similar dynamics. Both experiment and simulations show 
that the pit propagation speed and the precise shape of 
the knife edge ridges are universal, independent of ini- 
tial conditions, as predicted theoretically. These dynam- 




FIG. 1. Surface evolution of a magnesium alloy under uniform 
irradiation by a focused ion beam, after 5, 12, 19 minutes. 
The surface initially contained small holes that grew from 
scratches. The imaged region is 30/^m, viewed at 52°. 



ics can be understood by a theoretical analysis of the 
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FIG. 2. Diameter of pits in experiment (left) and simulations 
(right), as a function of time. Each marker represents a dif- 
ferent pit (expts) or initial condition (sims); experiments are 
accurate to ±0.15/im. Best- fit lines are dashed and legend 
indicates their slopes. The initial conditions for the simula- 
tions were: (cross) h(x,y) = — Qe~^ x +y \ (circle) h(x,y) = 
_ 10e -(* 2 +y 2 )/i y (triangle) h(x,y) = -^ e -^^ )lr{6 ^ with 
r(6) = l + 0.5sin(46>), tan0 = x/y. 




3. Numerical simulations of Eq. |T| with yield function 
Q, B = 1/100, at times t = 0, 11, 19. The initial 

-(* 2 + y 2 ) 

condition was h(x,y) = hge . 
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equations in which the knife-edge structure arises as a 
particular traveling wave solution with a large basin of 
attraction. Because there is only one such solution, the 
dynamics are relatively insensitive to the initial condi- 
tions and a pre-patterned surface will evolve to uniform 
knife-edge ridges with the same slopes and radii of cur- 
vature. If one can learn to control the location of the 
ridges by solving an inverse problem, one could poten- 
tially make a desired target pattern out of the knife-edge 
ridges. 
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FIG. 4. Top: slope h x through a horizontal cross-section 
y = 0, at the center of the simulations in Figure [3] at time 
£ = 11 (left), t — 15.5 (right). The ordinate axis indicates 
slopes bo = 2.3, 6* = 4.7 that act as dynamical attractors for 
the one-dimensional traveling wave equation. Bottom: maxi- 
mum slope \h x \ as a function of time. An initial transient pe- 
riod during which the narrow initial condition adjusts to the 
undercompressive shock bo, is followed by collision of shocks 
where the maximum slope jumps to 6* . 



Our experiments were performed on an ingot of a mag- 
nesium AZ91D alloy, with nominal alloying element con- 
tent 9 wt.% Al and 1 wt.% Zn; however, the experiments 
were confined to the alpha (aluminum-poor) phase of a 
two-phase mixture. The surfaces were uniformly irradi- 
ated using a FEI dual beam Focused Ion Beam (FIB) - 
Scanning Electron Microscope (SEM) delivering 30 keV 
Ga+ to the surface in a background pressure of 1.4x 10~ 6 
mbar at room temperature. The incident ion beam was 
parallel to the surface normal and the ion beam current 
was 3 nA. The beam was rastered in a boustrophedonical 
scan across a pre-defined region of the sample surface, 
during which time the beam would dwell at each discrete 
location for 0.1 fis and then move rapidly to an adjacent 
location. Separation between adjacent locations was set 
to nominally 50% overlap, which in this case meant a 
75 nm center-to-center spacing for a 150 nm diameter 
beam. The current profile within the beam is believed 
to be roughly Gaussian. The irradiation was interrupted 
periodically so that the irradiated surfaces could be ob- 
served using in-situ SEM from both normal incidence and 
tilted 52 degrees. 

The surface topography was initially irregular due to 
metallographic polishing scratches, and some of these ir- 
regularities initiated small holes in the surface (Figure [l] 
top.) Most of the holes then decayed to a flat surface, 
but certain holes developed into pits that continuously 
expanded. What is notable is the pits appear to be iden- 
tical: they expand at the same rate, and have sides of the 
same slope. We measured the diameters of five different 
pits as a function of time from a normal view of the sur- 
face. These changed at an average rate of 0.24/im/min, 
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with all rates lying within 0.02/im/min of the average, 
well within the resolution of our measurement of the di- 
ameter (Figure[2j left). We are not able to quantitatively 
compare the slopes, but qualitative examination of Fig- 
ure [l] middle, bottom suggests they are also very similar. 

When two pits collided they created very steep, sharp 
ridge-like structures. There are three examples of such 
ridges in Figure [I] (bottom) . These are notable because 
the length scale that is created by the collision is much 
smaller than any contained in the initial condition. Qual- 
itatively, the steep ridges have similar slopes in each of 
the three cases. 

We now turn to numerical simulation of these struc- 
tures, using a partial differential equation governing 
the evolution of the surface height h(x,y,t) on the 
macroscale: 

ht + R(b) + B \7 ■ ( -~J== Vk) = 0. (1) 

This equation is derived from the widely-used Sigmund 
theory of sputter erosiorP, by expanding the sputter 
integral for surfaces whose curvature is much smaller 
than the lat eral s cale over which an ion deposite its ki- 
netic energjEHSl. Here R(b) = RoVT+¥Y(b) is the 
the average velocity of erosion of the surface as a func- 
tion of its slope b = \Vh\ (or equivalently the angle of 
the incoming ion beam), obtained from the yield func- 
tion Y(b) by multiplying by a dimensional factor Rq 
and a geometrical factor. The fourth-order term with 
magnitude Bq is a function of the surface curvature 

K = V- ( v / 1 ^ b 2 VhJ , and models additional smoothing ef- 
fects such as Mullins-Herring surface diffusior l 16 l 17 lor ion- 
enhanced viscous flow confined to a thin surface layeaP^l 
We neglect t he se cond-order (curvature) terms that are 
often included^*', as the dynamics we are interested in 
occur when these are small. 

We have developed an efficient, stable method to solve 
Eq. (JlJ in two dimensions. Solving such higher-order 
nonlinear equations in multiple dimensions is a generally 
a challenge - explicit methods impose severe restrictions 
on the time step At for the scheme t o be s table (At < 
O(Ax) 4 , where Aa; is the grid spacing^212U) , while fully 
implicit schemes, that are unconditionally stable, require 
solving a difficult nonlinear problem at each time step. 
We overcome these difficulties by adding a fourth-order 
linear term to the equation that we treat implicitly, while 
treating the nonlinear parts of the equation explicitly^. 
Specifically, our scheme takes the form 

h J+1 +MAtA 2 h j+1 = 

h? - At (R(\Vh j \) - N{\Vh j \) + MA 2 h j ) (2) 

where h?(x,y) is the solution at time jAt, N(b) is the 
fourth-order nonlinear term, and M > is a con- 
stant that we are free to choose. Analysis for similar 
equations^HUl nas shown the scheme is stable irrespec- 
tive of the grid spacing provided M is large enough. How- 
ever, with fixed timestep the accuracy decreases if M is 



too large. We found a good balance between stability 
and accuracy with M = 1. 

To apply ([2]), the right-hand side is evaluated explic- 
itly using centered finite differences for the spatial deriva- 
tives, the result is converted to Fourier space using peri- 
odic boundary conditions, and /ij+i is found in Fourier 
space by solving the linear inversion problem. 

We perform simulations using a yield function of the 
Yamamura formal 

^gj = (cos^)" / ex P {-S((cos0 opt )~ 1 ~ 1)}, (3) 

where the parameters are / = 2.36, 9 opt — 69.5, and 
S = f cos 9 opt . Yamamura has shown that a great many 
experimentally measured yield functions can be repre- 
sented in this form, by fitting for /, 9 op t- Our theory (to 
be described) shows that the qualitative features of the 
dynamics are robust to changes in these parameters so 
we chose these for numerical convenience; we do not at- 
tempt quantitative comparison as the yield function for 
the experiments is unknown. 

We ran two kinds of simulations. First we looked at the 
formation of pits, by initiating the surface with individ- 
ual small holes. For small perturbations the pits expand 
slowly and decay without creating smaller length scales. 
However, for large enough perturbations, the pits evolve 
to a circular crater whose sides have a uniform, steep 
slope, that expands outwards with a constant velocity. 
The slopes of the sides, and the speed of propagation, 
are fixed numbers, independent of the form of the initial 
perturbation. Figure [2] (right) illustrates the constant 
speed for three different initial conditions. 

Next, to investigate how pits collide, we initiated the 
surface with two nearby pits. These expand, and when 
they collide, they form a ridge whose sides are very steep 
- much steeper than the sides of the crater. The slope 
of the sides is always the same, regardless of the initial 
conditions. Figure [3] shows the surface evolution for one 
choice of initial condition, and Figure [4] (bottom) shows 
the maximum slope as a function of time for this simu- 
lation. The plateau from time 13-20 corresponds to the 
slope of the ridge, and is the same height for a broad 
class of initial conditions. 

The simulations are striking because of their remark- 
able resemblance to the experiments. Two notable fea- 
tures occur in both: (1) Craters expand with a constant 
slope and velocity; (2) When craters collide, they create a 
very sharp knife-edge-like ridge, with steeper slopes than 
those originally on the surface. In both cases the slopes 
are universal, independent of the initial condition. 

Th ese fea tures were predicted by a recently-developed 
theorjUJMM ag ^ ne conse quence of the unusual type of 
traveling wave solutions that occur in the governing equa- 
tions (1). When a pit is large enough, the crater rim is 
locally nearly straight, and can be well approximated by 
a traveling wave that is invariant in one horizontal di- 
rection. Therefore we look for traveling wave solutions 
to the one-dimensional equation. As shown in Holmes- 
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Cerfon et. af^, the slope b = h x can propagate as a 
traveling wave provided the slopes in the far-field are 
held constant, so we look for solutions b = S(x — ct) to 

c{ S-b r )-(R(S)-R(b r)) = B 

(4) 

with boundary conditions 5 (+00) = b r , S(—oo) = —bi, 
S'(±oo) — 0. This equation is obtained from ([I]) by 
differentiating once, and then integrating from +oo to 
x — ct, while the speed c is determined by conservation 
of mass to be c = 1^1^ ■ Only certain pairs (bi, b r ) 
yield a solution, and this fact is crucial for determining 
the dynamics. 

Ahead of the crater rim the surface is flat, so we look 
for solutions with boundary condition b r = 0. The the- 
ory shows that there are solutions (6;,0) for all bi less 
than a critical slope b$ depending on the yield function. 
Above this critical slope there is exactly one boundary 
value yielding a solution: bi — bo, where bo is a number 
that again depends on the yield function. This solution 
is isolated and serves as an attractor for the dynamics in 
the following sense: if a one-dimensional surface is pat- 
terned initially to contain slopes greater than 6§, then 
these slopes will evolve spontaneously to the traveling 
wave connecting b to as the sloped region propagates 
into the flat far-fielcP^. 

This explains the first observation. When the ini- 
tial perturbation is large enough, it evolves to a crater 
whose rim propagates outwards with speed c correspond- 
ing to the discrete traveling wave (bo, 0), and whose sides 
therefore have slope bo. Figure [4] (left) shows the slope 
b through a horizontal cross-section of the simulated 
craters just before they collide, which clearly shows the 
slope of the sides is a constant, uniform value bo- 

What happens when pits collide? The simulations sug- 
gest a symmetry so we look for solutions with bound- 
ary conditions (bi,—bi). Again there is a threshhold de- 
termining the behaviour: when bi < b% for some yield- 
function-dependent number b% , there is always a solution, 
but when bi > b°, there is exactly one solution: bi = &*. 
This solution corresponds to a ridge with very steep sides 
and a small radius of curvature at the tip. Simulations^ 
showed that this solution will evolve from two nearby re- 
gions with slopes of opposite signs, provided both have 
(not necessarily equal) magnitudes > b%. Since bo > b%, 
we predict that colliding pits will evolve to the knife-edge 
ridge. 

Indeed, our numerical simulations confirm this - Fig- 
ure [4] (right) plots the slope b through a cross-section in 
the center of the colliding pits. The step-like appearance 
captures both of these traveling waves: the first step at 
bo is the slope of the original crater sides, and the second 
step at 6* is the knife-edge ridge. At later times (not 
shown) there is only one step, at &*, as the crater sides 
have entirely evolved to the knife-edge. 

For the yield function in our simulations the numerical 
values are b c = 1.26, b Q = 2.3, b% = 1.28, 6* = 4.7 - but 



the values of the dynamical attractors increase with opt , 
and for certain materials we predict ridges with slopes of 
b* = 30 or more^. Therefore we can create very sharp 
features by choosing an energy level or material that gives 
the desired values. 

' We have shown that the sharp, small-scale structures 
observed in our experiments and numerical simulations 
can be explained through the set of traveling wave solu- 
tions to the governing macroscopic equations. It is no- 
table that only the macroscopic mechanisms of erosion 
and smoothing are required to instigate the observed fea- 
tures. Of course, additional small-scale physics may help 
to explain some of the qualitative differences between the 
experiments and simulations, such as the curvatures of 
the pit bottoms. Indeed, we hypothesize that the experi- 
mental geometry may be significantly influenced here by 
multiple scattering effects, where ions incident on the pit 
wall and forward-scattered, as well as forward-sputtered 
atoms from the pit wall, may contribute to enhanced ero- 
sion along the perimeter of the pit bottom. 

The theory predicts that two traveling wave solutions, 
both with steep slopes, control the dynamics over a wide 
range of initial conditions. What is potentially useful 
about these solutions is that they arise spontaneously 
from smaller slopes - therefore we don't need to start 
with steep, small-scale structures in order to create them; 
these are created by the dynamics. This suggests a po- 
tential self-organizing principle for fabricating small-scale 
features on a surface, by pre-patterning the surface on 
the macroscale so that it evolves to a structure built of 
small-scale ridges. One is then interested in the inverse 
problem: to find an easily-achievable initial patterning of 
the surface, so that it evolves under uniform irradiation 
to a target small-scale pattern. 
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